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We study metal-insulator transitions between Mott insulators and metals. The transition mech- 
anism completely different from the original dynamical mean field theory (DMFT) emerges from a 
cluster extension of it. A consistent picture suggests that the quasiparticle weight Z remains nonzero 
through metals and suddenly jumps to zero at the transition, while the gap opens continuously in 
the insulators. This is in contrast with the original DMFT, where Z continuously vanishes but the 
gap opens discontinuously. The present results arising from electron differentiation in momentum 
space agree with recent puzzling bulk-sensitive experiments on CaV03 and SrV03. 
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Mechanisms and nature of correlation-induced metal- 
insulator transitions (MIT) are fundamental challenges 
in condensed matter physics*. Although the metals and 
insulators far away from the MIT are relatively well un- 
derstood, electronic states dominated and controlled by 
the proximity of the MIT are far from complete under- 
standing. We find many challenging phenomena such as 
the high-T c superconductivity in this region, which waits 
for better understanding of underlying proximity of MIT. 

The dynamical mean-field theory (DMFT) offers a 
clear picture of the MIT by taking the limit of high di- 
mensions 2 . In the DMFT, the MIT is approached from 
metals by the reduction of the quasiparticle weight Z at 
the Fermi level and the transition is driven by vanish- 
ing Z A . At the transition, a nonzero insulating gap is 
already open in the density of states (DOS) as a result 
of the vanishing Z of the coherent peak isolated from the 
well-separated upper and lower Hubbard bands. Early 
photoemission spectroscopy (PES) and inverse-PES of 
Cai-^Sr^VO^ showed a three peak structure, called a 
sharp coherent band at Ep and the incoherent ones cor- 
responding to the upper and lower Hubbard bands at a 
few electron volts above and below Ep, which is con- 
sistent with the above DMFT scenario. However, re- 
cent more bulk-sensitive PES have revealed a qualita- 
tively different and puzzling feature with a new broad 
peak at around 0.2eV connected to the lower Hubbard 
band around I.6eV and a pseudogap formation around 
the Fermi levels. This is in marked contrast with the 
DMFT results. 

In this paper, we show that a scenario of the MIT com- 
pletely different from the original DMFT emerges in two 
and three dimensional systems by extending the DMFT 
to allow the momentum dependence of the self-energy. 
The MIT is now not driven by a continuous reduction of 
Z to zero, but by the quasiparticle poles moving away 
from the Fermi surface, leading to a discontinuous jump 
of Z to zero at the Fermi level. Such a MIT is the conse- 
quence of an anisotropic extinction of the Fermi surface 
with the topological change, namely the electron differ- 
entiation in momentum space. We show the results by 
considering a cluster extension of the DMFT up to four 



sites for the Hubbard model on the square and cubic lat- 
tices. Our results indicate an opening of pseudogap in 
the DOS already in the metallic side leading to the con- 
tinuous opening of the true gap in the insulator at a finite 
energy resolution. This offers a picture consistent with 
the recent bulk-sensitive PES 5 -. 

In spite of the substantial success, limitations of the 
original (single-site) DMFT due to the ignorance of the 
momentum dependence in the self-energy have been re- 
peatedly pointed oui£i Recently, such as cellular DMFT 
(CDMFT)£, dynamical cluster approximation 9 , cluster 
perturbation theory^ -, self-energy functional approach 11 
and correlator projection theoryi^ have been developed 
with some promising results on the two-dimensional Hub- 
bard model by focusing on either doping or temperature 
effect£ ' 13 i 14 ' 15 i 16 ' 17 i 18 . Among these results, only Kyung, 
et. alJ£ first mentioned a four-band structure at half 
filling as the Mott transition is approached, which is 
remarkably different from the band structure predicted 
from the standard DMFT—. However, since no further 
analysis was performed, such as the analysis of the k- 
dependence of quasiparticle weight, self energy and spec- 
tral density, as well as the correlation functions, the 
mechanism and the criticality of the MIT cannot be de- 
rived from their results, while the new PES 5 casts doubts 
on the widely accepted MIT scenario^ and requires a new 
description of the MIT. Moreover, it remains unclear if 
the choice of cluster size and the dimensionality will affect 
the results. In this paper, we study in detail the two- and 
three-dimensional Hubbard model with the transfers for 
the nearest neighbor pairs with the standard notation: 

H = ~tJ2 (4r c jV + H - c -) + U J2 n iA n i,i- (!) 

on square and cubic lattices at half-filling with CDMFT. 
The CDMFT has aleady passed several tests against 
the exact results and the density matrix renormalization 
group method in one dimension^, where the CDMFT 
scheme is expected to be worst. We employ exact di- 
agonalization (ED) as an impurity solver of CDMFT 20 . 
The advantage of ED is the access to zero temperature 
and real frequency, as well as the ability to treat the 
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large interaction regime, which is hardly accessible by 
the quantum Monte Carlo method. 

Throughout this paper, we restrict our studies in the 
paramagnetic state. The magnetic or charge density 
wave order are indeed absent in SrVC>3 and CaVC>3. 
We take into account the temperature effect as well as 
the finite energy resolution by introducing a Lorentzian 
broadening factor. Finally, we should point out that, 
the CDMFT always underestimate the /c-dependent be- 
havior of self-energy owing to the limitation of cluster 
size. Therefore the possible non-Fermi-liquid behavior at 
asymptotically low-energy scale is hardly tractable. Nev- 
ertheless, even from our small cluster CDMFT study, we 
find that the scenario of MIT from DMFT breaks down. 

To implement the CDMFT approach, it is convenient 
to consider the cluster effective action derived from a 
cluster-impurity Hamiltonian of the form 

H ci = t »» c U c,/ ° + U 12 n ^ n »<i 

+ ^ ] £ma a l ncT a ma + ^ ] Vm^a( a \n<r c ti.o + H.C.), (2) 
ma m/jcr 

where the indices fi,v = 1,...,N C label sites within the 
cluster, while m = 1, ...,iV(, represent the bath degrees 
of freedom, and c M<T (a ma ) annihilates an electron on the 
cluster (bath), respectively. i M „ is the hopping matrix 
within the cluster, e m<J is the bath energy and V^„ CT is 
the bath-cluster hybridization matrix. We begin with an 
initial guess for the bath parameters e ma and deter- 
mining an initial cluster Weiss field Gg c (iuj). With this 
guess the cluster Green's function (GF) G° (iuj) is cal- 
culated by solving the cluster-impurity Hamiltonian ([2]). 
Applying the Dyson equation, = (Gq c )~ — , 
we obtain the cluster self-energy. To close the self- 
consistency loop, we obtain a new Gq' c (ilu) from 
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Here t Ik) is the hopping matrix for the superlattice with 



the wavevector k because of the inter-cluster hopping and 
D=2 for square lattice, 3 cubic lattice. The bath param- 
eters for the next iteration is determined by minimizing 
the following distance function 
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In the present study, we choose the size of cluster from 
1 up to 4. For 1-site CDMFT, we can reproduce the 
results of single-site DMFT. We set total system size N c + 
Nb > 10, i.e., we choose 10 bath sites for 1-site case and 
8 for all the rest cases. The self-consistency condition Q 




FIG. 1: (Color online) Density of state of Hubbard model 
on half-filled square lattice at four different values of U, ob- 
tained from CDMFT with cluster size up to 4. (a) is the 
1-site CDMFT which reproduces the results of the single-site 
DMFT. The central coherent peaks are narrowed but persist 
until the phase transition. (b),(c),(d) are 2-, 3-, and 4-site 
CDMFT, respectively. In all the three cases, the central peak 
is suppressed while a twin-peak structure appears around 
Fermi level. We take the broadening factor O.lt, viewed as 
a temperature effect or a finite energy resolution in experi- 
ments, a pseudogap is observed prior to the Mott transition. 



is imposed on the imaginary frequency axis with a cutoff 
frequency larger than all the bath energies. We set an 
effective inverse temperature (3t = 200 which serves as a 
low-frequency cutoff. The number of bath parameters is 
reduced by considering the symmetry of the square and 
cubic lattices. All the energies are measured in units of 
t. 

We first show in Fig. Q] the DOS for the half-filled 
square lattice. From Fig. Q] (a) to (d) , different choices of 
the cluster sizes up to 4 sites are presented. For single-site 
DMFT in Fig. [1] (a), the central coherent peak is pinned 
at the Fermi level and becomes narrow as Hubbard U 
increases. At U = lOt, it is shown that the central peak is 
well seperated from the upper and lower Hubbard bands. 
As U is larger than the transition point at U ~ lit, the 
quasiparticle spectral weight at Fermi level vanishes and 
the system changes from a Fermi liquid state directly to 
an insulator. The quasiparticle residue 



Z k = 1- 



duj 



(5) 



as function of increasing U is shown in Fig. [2] (a), where 
the lattice self-energy Y>i aU is equal to the cluster self- 
energy for the 1-site case. Therefore Z k is fc-independent. 
The quasiparticle residue continuously vanishes at the 
MIT which confirms the scenario of the MIT of DMFT-. 
However, as the cluster size become larger than 2, we will 
show later that the situation dramatically changes due to 
the involvement of nonlocal correlations. 
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FIG. 2: (Color online) (a) Momentum and U dependences 
of quasipaticle residue at Fermi surface for both cubic and 
square lattices obtained from 4-site CDMFT approach. The 
Lorentzian broadening factor is O.lt. The 1-site CDMFT re- 
sult for square lattice is also shown by the dashed curve for 
comparison, (b) The local spin correlation (open square), 
nearest-neighbor spin correlation (open circle), next-nearest- 
neighbor spin correlation (up triangle) and double occupation 
number (down triangle) as a function of U for square lattice 
obtained from 4-site CDMFT approach. Around (7 C =3.5 (5), 
the central peak splits into twin peaks for square (cubic) lat- 
tice, leading to a sudden jump of Z to vanish and a singular 
behavior in each correlation function. 



FigureQ](b) to (d) show that, while broad central peaks 
still exist for smaller U, they are suppressed as new twin 
peaks appear just below and above the Fermi level when 
U becomes larger (see U = 4i results calculated by 2-/4- 
site CDMFT and U = 6 At by 3-site one). Here we intro- 
duce a Lorentzian broadening factor t] — O.lt to mimic 
the finite energy resolution in experiments or the finite 
temperature effect. The observed formation of a pseudo- 
gap prior to the Mott transition suggested in the present 
result offers a very different scenario of the MIT from 
that of DMFT^. With further increasing value of (7, the 
new twin peaks are separated further and eventually a 
Mott gap smoothly opens (see U = 8t (7t) for 2-/3- (4-) 
site CDMFT respectively). 

The spectral density A (k,ui), the real ReE (k,ui) and 
imaginary ImE (k,ui) parts of the self-energy are shown 
in Fig. [3] (a), (b), and (c), respectively, at the Fermi 
surface X point k = (n, 0) for a 4-site cluster at various 
values of U on square lattice. The lattice self-energy E; att 
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At the X point, the slopes of ReE (fc, a;) are positive at 
U = At and 5t , where we observed a formation of pseudo- 
gap in the A(k,ui). Simultaneously, ImE (k,uj) show lo- 
cal minima at uj = 0, which indicates a non-Fermi-liquid 
behavior. These results signal the appearance of two new 
solutions in the quasiparticle equation uj — — ReE (fc, uj) 
in addition to the strongly damped solution at ui = 
which is also present in the noninteracting system. This 
is the consequence of the involvement of the short-range 
nonlocal correlations. We find that, with increasing U, 
not only the local magnetic moment is enhanced and the 
double occupation number decreases but also the nonlo- 
cal magnetic correlations within the cluster are enhanced 
gradually as shown in Fig. [5] (b). The nearest-neighbor 
spin-spin correlation is always negative so that the spin 
correlation is antiferromagnetic. It should be noted that 
the pseudogap appears even for the present cluster size 
where spin correlations are allowed to extend at most 
only up to the cluster size. 

The spectral density at the Fermi surface point fc = 
(7r/2,7r/2) (not shown) indicates that the pseudogap ap- 
pears more or less simultaneously for all the fc points 
at the Fermi surface for our limited cluster size up to 4 
sites. The more or less simultaneous pseudogap opening 
is presumably due to the smallness of the cluster which 
cannot allow for the nonlocal correlations larger than the 
cluster size. It is plausible that the pseudogap opens at 
different U for different fc points for a sufficiently large 
cluster. In the present results, we find that the separa- 
tion of the two new poles at fc = (tt/2, 7r/2) is always 
smaller than that at fc = (it, 0) as U increases. Such 
a differentiation among different fc points is already seen 
even in the Fermi-liquid region. In Fig. 13(d), we plot the 
spectral intensity A(k,ui) at w=0 in the Brillouin zone at 
U/t=3.3, Cold (hot) region is visible around (ir/2,ir/2) 
((7r,0)), implying a differentiation of MIT depending on 
the Fermi surface position. Near the MIT, the differen- 
tiation would make compatibly a gap in the hot region 
with a retained spectral weight in the cold region. For 
even larger U = 7t, a large separation of the new twin 
poles at each fc point with a strong peak in ImE (fc, uj) 
at the Fermi level leads to a large gap in A(Je,w). For 
smaller U — t and 3t, a Fermi liquid still survives. 

From Fig. [TJ we find that a gap opens at a smaller U /t 
for 2-, 3- and 4-site CDMFT than for the 1-site CDMFT. 
This at least partially originates from a cooperative role 
of short-range spin correlation for the gap widening. The 
Mott transition occurs at substantially lower U if the spa- 
tial antiferromagnetic correlation is taken into account. 
This indicates that the Mott mechanism of the gap open- 
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FIG. 3: (Color online) The spectral density A(k,u>) (a), real 
ReE (fc, w) (b), and imaginary ImE (k, uj) (c) parts of the self- 
energy at the Fermi surface X point k = (tt, 0) for a 4-site 
cluster at various values of U on square lattice. As U in- 
creases, a pseudogap develops in A(k,u>). A positive slope 
in ReE (k, uj) and a local minimum in ImE (k, uj) at u) = 
indicate the breakdown of the Fermi-liquid behavior. (d) Den- 
sity plot of A(k, 0) on square lattice in the metallic region at 
t/=3.3i.The Lorentzian broadening factor is O.lt 
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FIG. 4: (Color online) (a) and (b) are DOS of Hubbard model 
on half-filled cubic lattice, obtained from CDMFT approach 
with cluster size up to 4. (a) is the 1-site CDMFT which re- 
produce the results of single-site DMFT. (b) is 4-site CDMFT. 
The central peak splits into a twin-peak structure around 
Fermi level similarly to the square lattice. Here we take the 
broadening factor O.li. 



ing at the metal-insulator transition is cooperatively en- 
hanced by the Slater mechanism ascribed to the antifer- 
romagnetic order. 

As is seen from Fig. [TJ there still exists a small dif- 
ference between the results of 3 sites (odd number) and 
2-/4 sites (even number). The central peak persists to a 
larger value of U for the odd than for even sites. This 
may be due to the fact that the sharp coherent peak is as- 
sociated with residual entropy arising from the unpaired 
spin (spin doublet structure) induced in odd number of 
sites as is expected in 1- and 3-site clusters. 

In Fig. [2] (a) we also show the quasiparticle residue on 
the Fermi surface with increasing U obtained from the 
4-site CDMFT approach. The most prominent differ- 
ence from the single-site DMFT where the quasiparticle 
residue at the Fermi level vanishes continuously toward 
the transition point is that it decreases but remains at a 
nonzero value until new twin poles appear in the vicin- 
ity of the Fermi level at each k point. Then the MIT is 
characterized by escape of the quasiparticle poles away 
from the Fermi level resulting in a discontinuous jump 
of Z to zero, instead of the complete transfer of spectral 
weight from the well-seperated central peak to the upper 
and lower Hubbard bands. This is consistent with the 
recent study by the correlator projection theory^. From 



the CDMFT approach, the quasiparticle residue becomes 
/c-dependent along the Fermi surface. The residue in- 
creases monotonically from (ir, 0) to (tt/2, 7r/2) at a given 
U. Therefore, we only show two k points at Fermi surface 
in Fig. H](a). In Fig. [21(b), no jump is detected in all the 
correlations, indicating that the splitting of the quasipar- 
ticle poles at the Fermi surface into the twin poles in the 
vicinity of it is continuous. 

It is interesting to see whether the results qualitatively 
change if the dimensionality increases to three. We show 
in Fig. 2] the DOS for half-filled cubic lattice by using 1- 
site and 4-site CDMFT approaches. No qualitative dif- 
ference is found between square and cubic lattices. A 
twin-peak structure around the Fermi level is found by 4- 
site CDMFT approach as U increases (see U — 7t). The 
fc-dependent quasiparticle residue, as shown in Fig. [2] 
(a), remains nonzero at the proximity of the point where 
the Fermi-liquid behavior breaks down while the value is 
smaller than that of square lattice indicating a crossover 
to infinite dimension where the quasiparticle residue at 
Fermi level should continuously vanish. 

The results on the cubic lattice are relevant to the 
recent high-resolution bulk-sensitive PES of SrV03 and 
CaV03. At finite temperature T = 6K and with the en- 
ergy resolution 8meV, a peak around 0.2eV is observed 
with a pseudogap structure in qualitative but essential 



■5 



agreement with our results near the Mott transition. 

In summary, we have shown that the MIT scenario 
obtained by CDMFT with cluster size larger than 2 is 
quite different from that of the single-site DMFT. With 
the increase of the interaction U in the CDMFT, the 
quasiparticle peak becomes continuously split into twin 
peaks moving away from the Fermi level and leaving a 
pseudogap at the Fermi level on both the square and cu- 
bic lattices. With further increase of U, a full gap opens 
gradually from hot to cold spots. As a consequence, the 
quasiparticle weight remains fc-dependent but nonzero 
until the splitting of the central peak which leads to a 
sudden jump of Z to vanish. This is in contrast to the 
single-site DMFT, where it is /c-independent and vanishes 
continuously at the MIT point. We find that the transi- 
tion is not driven by a uniform and continuous vanishing 
of quasiparticle weight. The transition is rather driven 
by the escape of quasiparticle poles away from the Fermi 



surface. The momentum differentiation of the escaping 
process leaves a truncated structure of the Fermi sur- 
face such as the "arc" -type structure. This leads to a 
topological change of the Fermi surface at the transition 
to the insulator. The present unravelled mechanism of 
MIT characterized by the topological and momentum- 
resolved change in the Fermi surface indicates a route to 
the Mott transition without an isotropic divergence of 
the electronic effective mass and opens a way to study 
profound implications expected for the proximity from 
the MIT. 
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